Decoherence as a Probe of Coherent Quantum Dynamics 
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The effect of decoherence, induced by spontaneous emission, on the dynamics of cold atoms periodically 
kicked by an optical lattice is experimentally and theoretically studied. Ideally, the mean energy growth is 
essentially unaffected by weak decoherence, but the resonant momentum distributions are fundamentally altered. 
It is shown that experiments are inevitably sensitive to certain nontrivial features of these distributions, in a way 
that explains the puzzle of the observed enhancement of resonances by decoherence [Phys. Rev. Lett. 87, 074102 
(2001)]. This clarifies both the nature of the coherent evolution, and the way in which decoherence disrupts it. 



PACS numbers: 05.45.Mt, 03.65.Sq, 32.80.Lg, 42.50. Vk 

The theory of coherent quantum transport in periodic poten- 
tials is basic to solid state physics, and to our understanding of 
various conductance phenomena in crystal lattices. After suc- 
ceeding in isolating and manipulating single quantum objects 
such as ions or atoms, quantum opticians, as well as meso- 
and nanoscientists, have now started to build extended struc- 
tures of atoms or ions of increasing complexity. A natural 
way of doing so is to arrange one, two, or three dimensional 
regular aiTays of (cold or ultracold) atoms in optical lattices 
fl], which then can be considered as faithful realizations of 
strongly idealized, fundamental models of solid state theory. 
Beyond illustrating such theoretical models under clean and 
virtually perfectly controlled laboratory conditions, these ex- 
periments often also hold unexpected surprises (due to appar- 
ently innocent, real-life modifications of the original model), 
and promise highly rewarding technical applications in the fu- 
ture. Proposals that suggest using optical lattices for quantum 
information processing (2|] are but one example of this. 

In all such respects, the impact of noise and decoherence 
is a crucial issue fs', m', '5], because decoherence is expected 
to impair manifestly quantum phenomena. The present Re- 
port addresses a striking, seeming violation of this rule, which 
was experimentally observed with kicked cold atoms sub- 
jected to a pulsed, one-dimensional, spatially periodic opti- 
cal lattice la]. Here, a peculiar type of coherent quantum 
transport, called 'quantum resonance' is theoretically 

predicted for kicking periods rationally related to the propa- 
gation time of kicked atoms across the lattice constant. The 
mean kinetic energy of an atomic ensemble is then predicted 
to increase linearly with time, in sharp contrast to the be- 
havior predicted for non-resonant values of the kicking pe- 
riod, where it saturates (in the process of 'dynamical localiza- 
tion' 1 9], closely analogous to Anderson localization in one- 
dimensional disordered solids 1 10]). In previous experiments, 
enhanced transport was indeed observed for the lowest-order 



resonances fl fTHl . When decoherence was added to the ex- 
periment by controlled spontaneous emission (SE), the energy 
growth at resonance was found to be significantly faster than 
could be accounted for by the heating effect of momentum 
transfer due to SE i5 fl2ll . This looks like incoherent magnifi- 
cation of a purely coherent, and non-classical, phenomenon. 

In this Report we show how this counter-intuitive effect of 
decoherence can be resolved by inspection of the full atomic 
momentum distributions instead of merely their mean square 
(i.e., kinetic energy) values, as were considered in Ref. ]^. 
Decoherence then acts as expected: it destroys the coherent 
dynamics underlying the quantum resonances and, in particu- 
lar, certain nontrivial features of the momentum distributions 
^TJ. It is precisely this latter fact that produces the surpris- 
ing enhancement of the mean energy values observed in the 
experiment f^,'?']. 

Our experimental system is a realization of the 

paradigmatic Kicked Rotor (KR) model ]9, 13], extensively 
used in investigations of classical chaotic dynamics and its 
quantum counterpart fl^. After trapping and cooling in a 
magneto-optic trap (MOT), about 10^ cesium atoms are re- 
leased and, falling freely under gravity, are exposed to pulses 
from a vertical standing wave of off -resonant laser light. This 
is red-detuned from the 6^81/2 ^ B^Pi/j, [F = A ^ F' = 
3) DI transition by (5l = 27r x 30 GHz, and has a wavelength 
Xl = 894. 7 nm. On release, the atomic temperature is 5/iK, 
corresponding to a Gaussian momentum distribution with full 
width at half maximum (FWHM) 1 2 where A:^ = 27r/AL. 
The duration of each (square) pulse is tj, = 500 ns, and the 
peak intensity in the standing wave is ~ 5 x 10-* mW/cm^. 
Due to the ac Stark shift, these pulses result in (5-function- 
like applications of a sinusoidal potential, with spatial period 
Al/2. Classically, the maximum impulse that this can im- 
part is hG(t>d, where (t>d = ^l^tp/86L, and fl is the Rabi fre- 
quency of the atoms at the intensity maxima of the light field. 
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Quantum mechanically, it imparts momentum to the atoms 
in integer multiples of HG, where G — ik^. Both the den- 
sity distribution of the trapped atoms and the standing light 
wave intensity profile are Gaussian, each with FWHM 1 mm, 
so the mean value of (f)d as experienced by the atomic ensem- 
ble is ~ O.Stt. The standing wave passes through a voltage- 
controlled crystal phase modulator which can shift the posi- 
tion of the standing wave between consecutive pulses so that 
it effectively cancels the effect of gravity. Thus in the rest 
frame of the atomic ensemble the standing wave appears to 
be stationary, yielding KR dynamics with kicking period T 
(S> tp), despite the presence of gravity. Decoherence can 
be introduced by inducing SE in the atoms through applica- 
tion of an additional 2/xs pulse of laser light after each of the 
kicking pulses. This light is 60 MHz red-detuned from the 
6^Si/2 62p3/2, {F = A ^ F" ^ 5) D2 ti-ansition. The in- 
tensity of the pulse can be controlled so that the mean number 
of SE per atom per pulsing cycle, TTse, can be varied contin- 
uously. Finally, after application of the pulses, the atoms fall 
through a sheet of light resonant with the D2 transition, lo- 
cated 0.5m below the point of release. This allows us to deter- 
mine their momentum distribution by a time-of-flight (TOE) 
technique, with a resolution of 2± KG. 

In the absence of SE, the Hamiltonian that generates the 
time evolution of the atomic wave function may be written in 
the following dimensionless form 



(j>d cos(2;) S{t — mr) , 



(1) 



where p is momentum in units of KG, x is position in units 
of G^^, and M is the mass of the atoms. The units of energy 
and time are then li^G^ /M and A// HG^, respectively. In such 
units the kicking period is r = hG^T/M. This Hamiltonian 
is very close to that of the well-known (5-kicked rotor, with 
the sole (but important, as discussed below) difference that 
our cold atoms are moving along a line rather than in a circle. 
Notwithstanding this, the dynamics of the atoms reflect char- 
acteristic properties of the quantum KR (whose correspond- 
ing classical phase space is mixed, consisting of regular and 
chaotic components, for non-vanishing kicking strength lITsll V 
The nature of the quantum transport sensitively depends on 
the parameter r. If r = Anr/ q, with r, q integers, then the 
kicking period is rationally related to the propagation time of 
kicked atoms across the lattice constant, and quantum trans- 
port is typically enhanced by quantum resonance 1 8 ] . If t/Att 
is sufficiently irrational, then transport is inhibited by quan- 
tum interference, i.e., by dynamical localization |9]. 

A conceptually simple way of experimentally testing this 
theoretical picture is to measure as a function of r the mean 
kinetic energy of the atoms, henceforth to be referred to as 
'mean energy', after a fixed interaction time of N kicks. The 
result of such a measurement at = 30 in the absence of 
SE is shown in Fig. [Ha), with O.lOvr < r < G.SItt. The 
mean energies were extracted from a finite momentum win- 
dow (—60 < p < 60), and low-amplitude noise in the time- 
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FIG. 1: Experimental values of the mean energy of the atomic en- 
semble after 30 kicks, as r is varied from O.IOtt to 6.3l7r (i.e., 
6.5 fJ,s < T < 210.5/is) in (a) the absence, and (b) the presence 
of induced spontaneous emission, with Use — 0.14. 



of-flight signal was eliminated by imposing a signal threshold 
estimated from the background noise level at high momenta. 
Figure^a) shows an underlying smooth, periodic dependence 
on r. This reflects the r-dependence of the localization length 
of the dynamically localized atomic sample 116]. The struc- 
ture superimposed on the basic periodic variation has partic- 
ularly narrow peaks at r = 2tt, Att, 6tt. These are the main 
quantum resonances, with q = 1, r = 1 and g = 2, r = 1, 3. 
Higher-order (q > 3) resonances were not unambiguously re- 
solved within the given observation time. The present Re- 
port is therefore focused on the main resonances, and specifi- 
cally at the surprising way in which they react to decoherence, 
as shown in Fig. ^h). The mean energy growth at the reso- 
nances is clearly enhanced. In other words, resonant transport, 
which is due to constructive quantum interference, appears to 
be stabilized by decoherence (while the dynamical localiza- 
tion away from resonance is barely affected, confirming that 
we are in the regime of weak decoherence |3]). The apparent 
inconsistency of these experimental observations with what 
seems theoretically reasonable is resolved as follows. 

An atom periodically kicked in space and time is described 
by a wave packet tp{x) composed of 27r-periodic Bloch states 
-0/3 (x), that is. 



dp exp(i/3a;)V'/3(a;) . 



(2) 



where f3 is the quasimomentum. In our units, it is given by 
the fractional part of the momentum p — n + f3 (n € Z). It is 
conserved in time, so the different Bloch states in (|2j evolve 
independently of one another, and their momenta only change 
by integers. Under the resonance condition r — Airr/q, a spe- 
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cial situation occurs for a specific, discrete subclass of values 
of p. Besides being periodic in coordinate space, the one- 
period evolution (Floquet) operator is then also periodic in 
momentum space, with the integer period q. This happens 
when P = 'm/2r, < m < 2r, m integer. The amplitudes of 
such waves at momentum states separated by qhG (in physical 
units) exactly rephase after each kick |6]; here we specialize 
to q = 1,2. This rephasing is analogous to the Talbot effect 
in optics, so we speak of these resonances as occurring at ra- 
tional multiples of the half-Talbot time T1/2 = 2t: M/hG'^ = 
66.7/xs (for which r — 2tt). In much the same way as spa- 
tial periodicity enforces ballistic motion in physical space, the 
momentum periodicity which holds for special values of /3 
(i.e., /3 1/2 for q = 2, and = 0, 1/2 for g = 1) enforces 
ballistic propagation of the corresponding states in momentum 
space; thus their energy grows quadratically in time. The re- 
maining Bloch components of the original wavepacket, with 
/3 not in the 'resonant' class, undergo a quasiperiodic energy 
exchange with the driving field, leading to a finite spread of 
the associated (/3-dependent) momentum distribution for all 
times. Upon incoherently averaging over the continuous set of 
quasimomenta which constitute the atomic ensemble, there is 
competition between quasi-periodic and ballistic propagation, 
and as N increases the values of /3 that populate the ballis- 
tic growth must match more closely the ideal resonant values. 
On the one hand, this leads to linear growth of the total mean 
energy, E k, (jP^N/A. On the other hand, a stationary momen- 
tum distribution Ps{n) ~ limjv^oo Pin, N) (7|, given by 



da 
2^ 



'if it a)), (3) 



emerges, where N) is the coarse-grained momentum 
distribution of the ensemble after N kicks. The coarse- 
graining is on the scale of unity {HG in physical units) so 
as to yield a distribution in n, which is consistent with the 
finite-size binning of the experimentally detected momen- 
tum distribution. In Eq. Q, hin') is the initial (assumed 
smooth) momentum distribution, /(^, a) = (pd sin(^) csc(q;), 
^ = 7r(2/3 — 1), and Jn-n' is a Bessel function of first kind 
and order n — n'. The asymptotic distribution Ps (n), shown in 
Fig. 13 a), is attained because the phases a of the nonresonant 
Bloch components of the original wavepacket, accumulated 
under the action of the time evolution operator, are effectively 
averaged. 

For finite times, P(ri, N) exhibits a narrow, stationary peak 
centered around n ~ Q, algebraic decay oc over inter- 
mediate momenta, and 'ballistic wings' due to the almost- 
resonant (3 values, which move to higher momenta linearly in 
time. It is important to note that the central peak is narrower 
than the exponential distribution observed in dynamically lo- 
calized atomic ensembles, and that the linear energy growth is 
observed at all times, in spite of the onset of the stationary dis- 
tribution. There is no inconsistency here: the asymptotic limit 
to which the distribution tends has an fall-off at large n, 
and hence has an infinite mean energy (i.e., a divergent 2nd 
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FIG. 2: Normalized experimental (solid line) and coarse-grained nu- 
merical (dotted) momentum distributions after = 30 kicks at the 
quantum resonance r ~ 27r (T = 66.5/^s) in (a) the absence, and 
(b) the presence of induced spontaneous emission. Experimentally, 
risE = (0.14 ± 0.04), while tTse = 0.1 numerically. The dashed 
curve in (a) is the asymptotic distribution Ps{n), as given in Eq. js}- 
The arrow labels in (a) indicate the ballistic wings, whose momen- 
tum varies linearly with TV. Note the slight asymmetry in the exper- 
imental distributions around n = 0, due to non-ideal aspects of the 
realization. The fainter dotted lines show the signal threshold and 
momentum cuts imposed on the experimental data when calculating 
mean energies. 



moment). The ballistic wings, also experimentally observed 
in Ref. |lll], dominate the theoretically computed mean en- 
ergy growth. As the wings are fed by the resonant-/? subclass, 
the resonant energy growth is ultimately due to conservation 
of the quasimomentum (3 in the kicking process. 

Experimental detection of these wings in the final atomic 
momentum distribution is extremely difficult, for several rea- 
sons. The most important of these for our present discussion 
is that the wings must not have moved beyond the cutoffs im- 
posed by the signal-to-noise ratio. Though relatively small 
in terms of population, this experimental loss from the ballis- 
tic wings leads to a mean energy which is significantly less 
than the theoretical value. Our theoretical picture is compared 
with experimental data in Fig. |2ja), where experimental and 
numerical momentum distributions are shown at the quantum 
resonance r 2± 27r (T = 66.5/is), after = 30 kicks. Note 
that our numerical simulation exhibits the ballistic wings of 
the distribution, which are swamped by the noise background 
in the experimental data. When processing such data, only 
momenta in the window [—60, 60] were taken into account. 
Furthermore, the experimental distribution exhibits an asym- 
metry around n = which is not present in the theory. This 
is due to two effects: the first, and most important, is that of 
the lock-in amplifier and its associated low pass filter, used in 



4 



the TOF measurement, which slightly distort the momentum 
distribution. The second is that the removal of gravity's effect 
by the crystal phase modulator is imperfect; gravity breaks 
the symmetry of the system's evolution and hence of the mo- 
mentum distribution. Nevertheless, the experimental and the- 
oretical results agree very well in the central part of the dis- 
tribution. Other deviations of the experimental system from 
the ideal are: (i) pulses are not (5-like as they have a finite du- 
ration tp, (ii) random amplitude noise is introduced by laser 
power fluctuations (±5%), and (iii) different atoms are sub- 
ject to somewhat different values of 4>d 1 12]. 

The addition of noise reshuffles the quasimomenta of the 
initial distribution, at a rate proportional to tTse, and thus de- 
stroys the conservation of quasimomentum. This reshuffling 
prevents atoms from remaining in the fast-travelling quasi- 
momentum range for a long time; the formation of ballistic 
wings is thus inhibited. On the other hand, reshuffling gives 
all atoms a chance of sojourning a while in those quasimo- 
mentum ranges and hence experiencing a transient ballistic 
momentum growth. As a result, the distribution at moder- 
ate momenta broadens in time, at the expense of the ballistic 
wings. This is seen in Fig.|2jb): for tTse — 0.14, the momen- 
tum distribution is strongly broadened as it exhibits enhanced 
population of moderate momentum states. 

The incoherent dynamics are amenable to analytical treat- 
ment ||7i], which shows that the distribution no longer ap- 
proaches a stationary form. Instead, it evolves towards 
a continually broadening, diffusion-like Gaussian. The 
theoretically-obtained line shape compares very favorably to 
the experimental one, as shown in Fig. |3b). Note that we 
chose Use = 0.1 for the analysis in order to achieve an op- 
timal fit to the experimental data. The uncertainty in the in- 
tensity of D2 light experienced by the atoms, due to loss at 
the glass faces of the vacuum system and the exact shape of 
the 2fis pulse, means that the experimental value of Use is 
(0.14 ± 0.04), consistent with this best theoretical fit. 

Analytically, the mean energy grows according to £' ~ 
{D/2 + (t>1l^)N, where D ~ Use/I^ is the diffusion coef- 
ficient associated with the momentum transfer due to SE 0]. 
Since D is rather small (~ 0.01) for the cases considered here, 
the mean energy growth is almost the same as in the resonant 
case without decoherence, where the same expression for E 
is obtained, except D = (see above). Hence weak decoher- 
ence destroys the conservation of quasimomentum, which lies 
at the very root of resonances, yet in the ideal model it only 
affects the resonant energy growth mildly. However, the theo- 
retically almost identical energy growth is produced by quite 
different physical mechanisms, which react to the cutoffs in- 
herent in experimental detection schemes in dramatically dif- 
ferent ways. In the coherently-evolved case, the growth of 
the experimentally measured energy is strongly depressed as 
soon as the ballistic wings escape the detection windows. In 
the presence of noise, the energy growth is not due to the bal- 
listic wings, but rather to diffusive broadening of the whole 
distribution, and is therefore dominated by the center Hence 
the effect of finite experimental detection windows is much 



less severe and the measured mean energy remains closer to its 
ideal value up to a higher value of N. This leads to an apparent 
enhancement of the resonance peaks compared to the SE-free 
case. Such noise-induced signal enhancement is reminiscent 
of 'stochastic resonance', where the response of a system to 
some input signal is enhanced by stochastic activation |17|. 
However, the hallmark of stochastic resonance is a maximum 
signal enhancement at an optimal, nonvanishing noise level. 
This has not so far been established in our present scenario. 

In summary, we have shown that the linear growth with 
time of the mean energy at quantum resonance, inhibited in 
experiments by finite detection windows on finite time scales, 
is restored there by adding noise. This effect is ultimately 
rooted in the difference between atoms and rotors and is a 
striking, albeit indirect, demonstration of the peculiar nature 
of coherent resonant transport for kicked atoms, and how it is 
modified by photon recoil-induced decoherence. 
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